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Influence of sequence correlations on the adsorption of random copolymers onto 

homogeneous planar surfaces 

Alexey PolotskyQ Friederike Schmid, and Andreas Degenhard 

Fakultdt fiir Physik, Universitat Bielefeld, Universitatsstrafle 25, D-33615 Bielefeld, Germany 

Using a reference system approach, we develop an analytical theory for the adsorption of ran- 
dom heteropolymers with exponentially decaying and/or oscillating sequence correlations on planar 
homogeneous surfaces. We obtain a simple equation for the adsorption-desorption transition line. 
This result as well as the validity of the reference system approach is tested by a comparison with 
numerical lattice calculations. 
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I. INTRODUCTION 

Random heteropolymers (RHPs) currently attract 
much attention and have been studied intensely in the 
last two decades. The reason for the growing interest in 
them is twofold. On the one hand, they can be used to 
create new polymeric materials, or to modify and improve 
surface and interface properties. On the other hand, 
RHPs serve as simple model systems for biologically rele- 
vant heterogeneous macromolecules, such as proteins and 
nucleic acids. Studying the physics of RHPs helps to un- 
derstand generic properties of these macromolecules and 
phenomena which play an important role in living organ- 
isms, such as protein folding, denaturation/renaturation 
of DNA, secondary structure formation by RNA, etc. 

In the present paper, we consider RHP adsorption on 
homogeneous and impenetrable (solid) surfaces. This has 
technological relevance in the context of adhesion (glue) 
and surface design. Studying RHP adsorption can also 
be regarded as a (very) first step towards understanding 
some elementary principles of molecular recognition. 

For entirely random RHPs, the problem was investi- 
gated theoretically by a number of authors. Joanny 01 
considered the adsorption of polyampholyte chains on a 
planar solid surface. In the case that the electrostatic in- 
teractions are strongly screened and thus short-ranged, 
he treated the polyampholyte chain as an ideal RHP 
with no intersegment interactions, and solved the prob- 
lem using the replica trick and the Hartree approxima- 
tion. A replica symmetric mean field theory that takes 
into account intersegment interactions has been devel- 
oped by Gutman and Chakraborty 0j I3j to describe 
RHPs interacting with one or two parallel Q planar 
surfaces. They analyzed the interfacial structure (as ex- 
pressed through the density and the order parameter pro- 
files) as well as the adsorption-desorption phase diagram. 
Stepanow and Chudnovsky Q combined the replica trick 
with the Greens function technique to study the localiza- 
tion of a RHP onto a homogeneous surface, and suggested 
a phase diagram for the localization-dclocalization tran- 
sition. The system has been also studied by Monte-Carlo 
simulations 0, 0] ■ Sumithra and Baumgaertner |^ have 
complemented their simulations by a scaling analysis of 
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the conformational and thermodynamic characteristics of 
the system. 

In all of these works the RHPs were taken to be of 
BernouUian type, i. e., the monomer sequences had no 
correlations. Sequence correlations should, of course, af- 
fect the adsorption properties of RHPs. van Lent and 
Scheutjens have formulated a self-consistent field the- 
ory for random copolymers with correlated sequences in 
the case of annealed disorder, i. e., the copolymer se- 
quences can "adjust" to their local environment. They 
used it to study the adsorption of random copolymers 
from solution within a lattice model. Zheligovskaya et 
al. have explored the idea of conformation-dependent 
sequence design and compared the adsorption proper- 
ties of AB-copolymers with special "adsorption-tuned" 
primary structures (adsorption-tuned copolymers, ATC) 
with those of truly random copolymers. Monte-Carlo 
simulations revealed that some specific features of the 
ATC primary structure promote the adsorption of ATC 
chains, coinpared to random chains. Chakraborty and 
coworkers l^, ^3 have demonstrated the importance of 
sequence correlations in a number of studies on recogni- 
tion between random copolymers and chemically hetero- 
geneous surfaces. 

Whereas these studies were mostly numerical, De- 
nesyuk and Erukhimovich jl^ have recently presented an 
approximate analytical solution of a closely related prob- 
lem, the adsorption of RHPs with sequence correlations 
on penetrable homogeneous substrates (e. g., interfaces 
between two selective solvents). Their study is based on 
a reference system approach originally due to Chen [T5l |. 
The sequence correlations were described by a correla- 
tor introduced by de Gennes that covers both expo- 
nentially decaying and oscillating correlations. Denesyuk 
an Erukhimovich showed that chains with both types 
of correlation exhibit a continuous adsorption transition, 
and constructed the corresponding adsorption-desorption 
phase diagrams. 

In the present paper we extend this work to study the 
adsorption of RHPs with correlated monomer sequences 
onto a homogeneous impenetra ble p lanar surface. We use 
the reference system approach and follow closely the 
ideas developed in Ref. ui!. This will allow us to derive an 
explicit and surprisingly simple analytical expression for 
the adsorption transition point. To test our prediction, 
we also performed numerical calculations on the lattice, 
and compared the results to the premises and the result 
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of the analytical theory. 

The rest of the paper is organized as follows: The an- 
alytical theory is developed in Sec. ^ We begin with 
defining the model in Sec. Ill Al and discuss the theo- 
retical approach in the next three subsections, Sec. Ill Bl 
Sec. Ill CI and Sec. Ill Dl The final result is presented in 
Sec. Ill El Then, we describe the numerical calculations 
in Sec. IIIII The numerical calculations are compared to 
the analytical prediction in Sec. IIVI We summarize and 
conclude in Sec. 

II. ANALYTICAL APPROACH 

A. Definition of the Model 

We consider an ensemble of single heteropolymcr Gaus- 
sian chains, each consisting of monomer units, near 
an impenetrable planar surface. Monomers only interact 
with the surface, not with each other. The strength and 
the sign of the monomer-surface interaction varies from 
monomer to monomer and is characterized by an affin- 
ity parameter ^. A particular heteropolymcr realization 
is described by the sequence ^(n), where n indicates the 
monomer number (n = 1, 2, N). It is, however, conve- 
nient to consider n as a continuous rather than a discrete 



variable. The Hamiltonian of the system can then be 
represented as follows: 

where the z direction is perpendicular to the surface 
and we have integrated out the x and y coordinates. 
Here f3 = 1/kT is the Boltzmann factor as usual, a is 
the Kuhn segment length of the chain, and V{z) gives 
the shape of the monomer-surface interaction potential 
(/ V{z) dz = 1), which is taken to be the same for all 
monomers. More specifically, we use a pseudopotential 
- a Dirac well shifted by the small but finite distance zq 
from the impenetrable surface 

V(A - / '^(^- ^o) iiz>0 , . 

" \ +00 if z < 

The parameter zq > corresponds to the range of the 
surface potential. It has been introduced mainly for tech- 
nical reasons, which will become apparent further below. 
In the end (Sec. Ill E|l . we will take the limit zq — > 0. 

The sequences ^{n) are randomly distributed according 
to a Gaussian probability function 
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Here ^0 = (0 the mean affinity and c(ni,n2) describes 
the sequence correlations. The latter is taken to have the 
special form suggested by de Gennes 0| , 

c(ni,,i2) = ((e(^i)-Co)(e("2)-^o)> 

= c(|n2-rii|) = A2Rc e-^l"^-"il, (4) 

where Re denotes the real part of a complex number. This 
allows us, by proper choice of the parameter F, to cover 
both exponentially decaying (real F) or oscillating (purely 
imaginary F) correlations. The expression J I?{^(n)} de- 
notes functional integration over all possible sequences 
(integration in sequence space), and c~^(n) is the inverse 
function of c(n), defined by / dn' c~^{n — n')c{n') = S{n). 

Wc will calculate the quenched average over all se- 
quences, i. c, the free energy of the system is the disorder 
average of the logarithm of the conformational statisti- 
cal sum, F = (loffZ,„„f)ci-„\. For infinitelv long chains, 



the adsorption transition in this model is encountered at 
F 0. 

B. The Reference System Approach 

Having defined the model, we now turn to describing 
the analytical approach. To perform the quenched aver- 
age we exploit the well-known replica trick 

Z"^ — 1 

(log^)c(n) = lim ( )?(„)■ (5) 

This reduces the problem of evaluating the disorder aver- 
age over a logarithm, (logZ), to the much more feasible 
task of evaluating the average partition function for a set 
of m identical heteropolymcr chains (m replicas) (Z™). 
Inserting Eqs. (El - 101) and using properties of the Gaus- 
sian distribution, we obtain 
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where a is the repUca index and / D{z(n)} denotes the 
integration over aU possible trajectories z{n) of the chain. 

One can easily see that the problematic term in this ex- 
pression is the last factor which couples different replicas. 
To simplify this term we introduce a reference homopoly- 
mer system whose conformational properties are rea- 
sonably close to those of the original heteropolymer. It is 
natural to choose as reference system a single homopoly- 
mer chain with uniform monomer affinity S,o + Vq, which is 
exposed to a surface potential V{z) of the same shape (j2l 
as the original RHP. The parameter vq is a fit parameter 
which can be adjusted such that the reference system is 
as close as possible to the original system. It is nonzero 
and positive, because the effective surface affinity of het- 
eropolymers is generally larger than that of homopoly- 
mers with the same charge. For example, even globally 
neutral random heteropolymers with = can adsorb 



onto a surface, whereas neutral homopolymers obviously 
remain desorbed. 

The Hamiltonian of the reference system has the fol- 
lowing form: 



Ho 



kT 2a2 



TV 



dn 



dz 



dn 
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-mo + vo) dnV{z{n)) 
Jo 

(7) 

The difference between the free energies of the original 
and the reference system is given by 



-P{F-Fa) = lim 



- Z, 







(8) 



where is the partition function of the m times repli- 
cated reference system: 



J]l?{z„(n)}exp "^"(^j +/3(eo+«o)y^ dn V{z^{n)) 



(9) 



Thus one has 



n^^KNIxexp -^Y. dn(^) +P{^o+vo) dnY, V{z^{n)) 
exp — / dni dn2 y(za(ni)) y(z^(n2))c(7i2 - rii) - /3uo / dn V{za{n)) 

. ^ -^0 -^0 a=l 7=1 a=l 



- f 



(10) 



r 



The first exponential in this expression corresponds 
to the product of the m independent propagators 
'P{^0:i'o, Za{n)) for the m chains in the rephcated ref- 
erence system, which give the statistical weight of the 
particular trajectories Za{n). Since the original and the 
reference systems are assumed to be close to each other, 
we expand the exponent in the second term: — 1 ~ A. 



Here A is the difference between Hamiltonians of repli- 
cated original and the reference systems —(3{H™ — H^). 
This means that we approximate the difference between 
the free energies of the original and the reference system 
AF = Fq — F hy its first cumulant AFi . The expansion 
gives 



^"-^o"- J n [nzc.{n)}V{^o,vo,Za{n)) 



2 rN 



— / dni I (in2 F(2:a(ni)) V^(z-y(n2))c(n2 - ni) - /Jwo / dn^^ V{za{n)) 
Jo Jo ^ 1 Jo 1 

^ ^ a— 1 7=1 a— 1 



(11) 
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The second integral in Eq. Hll() is equal to 

/2=y n [T^{za{n)}r{^o,Vo,Zc.in)) 



a = l 

•^0 a=l 



(12) 



= PvoNmViz{n))Z^, 

where denotes the average with respect to the refer 
ence system: 

JV{zin)}Gi^o,vo,zin))X{z{n)) 



^ JV{zin)}Gi^o,vo,zin)) 
Next we consider the first integral 

/i = y n [T^{^^c.{n)}V{^o,vo,z^{n)) 



(13) 



Q = l 

2 



0- 



/ drii / (jr7,2 V{za{ni)) 
Jo Jo 



(14) 



X ^ ^^(2;7(n2))c(n2 - ni). 

7=1 

The product of the two sums in Eq. H14() gives 

m m 

^n^„(ni))^l/(z^(n2)) 

a— 1 7—1 

= Y,VMni))V{z^{n2)) 

m m 



(15) 



.^^F(z„(ni))F(z7(n2)), 



where we have separated contributions from the same and 
different replicas. This leads to the following result for 

h = f3^m ^o" / dk c{k) {N - k) (16) 



\{m-l)V{z^{n))V{z^{n + k)) 
+y(z„(n))y(z„(n + fc))}. 



a/7 



Finally, collecting all terms, taking the limit m ^ 0, 
dividing by » 1, and using c[k)k/N :s> 1 for all k, 
we obtain for the first cumulant AFi 



r/ dkcik)Viz^{n))V{za,{n + k)) 
-'^ Jo 



N 



/ dkc{k)V{zain))V{z^{n + k)) 
Jo 



(17) 



The next step is to find the Greens function of the refer- 
ence system, which is needed for calculating the averages 
in Eq. 



C. The Reference System 

As defined above, the reference system is a homopoly- 
mer chain with the affinity parameter vjq = ^o + ^Oi which 
propagates in the half space z > subject to the poten- 
tial V{z) = 5{z — zq). The Greens function of the poly- 
mer chain in the reference system satisfies the differential 
equation 



dG 
'dN 



a" d'^G 



+ l3wo5{z^zo)G{z,z'-N) 



6 dz'^ 

with the initial condition 

G{z,z']0) = 5{z- z') 

and the boundary conditions 

G(z,z';7V)U^oo = 
G{z,z';N%^o = 



(18) 



(19) 



(20) 



Eq. (|18|l has the structure of a Schrodinger equation and 
can be solved analogously by a separation of variables. 
The solution is an expansion in eigenfunctions of the 
"Hamilton operator" , 



G(z,z';^)=5]e--^^,(z)V'*(^'), 



(21) 



where Si and '4'i{z) are eigenvalues and eigenfunctions of 
the stationary equation 



a' ay 

6 dz"^ 



+ (3woS{z - zq)-\\){z) = -eV(-z)- 



(22) 



For very long chains (N ^ 1), the sum H21|) is domi- 
nated by the largest term, corresponding to the lowest 
eigenvalue (— eo) (ground-state dominance): 



G(z,z';iV)^e-«^7^o(^)^o(^')- 



(23) 



The analogous problem of quantum-mechanical parti- 
cle motion in the potential Viz) has been treated in detail 
by Aslangul |23| . To establish the connection to these cal- 
culations, we rewrite the stationary differential equation 
in the form 



ii"{z) - k'^ijj^z) = -koS{z - zo)')Jj{z), 
where k and fco are given by 

6eo 



and 



k' = -- 



ePwQ 6/3(^0 + Wo) 



(24) 



(25) 



(26) 



The boundary conditions are: Tp{z)\z=o ~ and 
^(z)U^ po = 0. Eq. has the same form as Eq. (3) 
in Ref. |2j except that we have chosen ko > 0, and we 
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can simply apply the results of Ref. |2l|. For the ground 
state eigenfunction we obtain 



g-fcb|z-2o| _ g-fcb(z + 2o) 



(27) 



where fc^, determining the ground-state energy 1)25(1. is 
the root of the transcendental equation 



2k 



-2k Zr 



= fcn 



(28) 



Since different replicas are independent, the second 
term of Eq. H17|) becomes 



V{z^{n))V{z^{n + k)) =V{zo,{n))' 



2fc2/fco 



1 + {2kb - ka)zo 



2 (33) 



and the normalization constant A is equal to 



A- 



fco/2 



1 + (2fcf, - ka)zo 



1/2 



(29) 



Eq. (|28|l always has the trivial solution fc = 0. A non- 
trivial solution exists if the slope of the line y = 2fc is 
smaller than the slope of the curve y = fco(l — e"^*^^") at 
z ~ 0. This condition yields the value of the adsorption 
transition point. 



l.tr 



1 

20 ' 



(30) 



The expression for the ground-state Greens function is 
fco/2 



G{z,z']N) = 



1 + {2kb - ko)zo 



-kh\z'-zo\ 



'k^{z'+za) 



(31) 



Based on these results we can now calculate the aver- 
ages in (|17ll . using the ground-state dominance approxi- 
mation 123() for the Greens function. 



V{zin)) (32) 
/cT Io° "^"^ ^(^' l)S{z — zq)G{z, y; N ~ I) 



dl J^^ dx dy dz G{x, z; l)G{z, y; N — I) 
dz\Mz)\^6{z~zo)^\Mzo)\^ 



= A^ (1 - e-^'^'^'Y = A^ 



2h\^ _ 2fc2/fco 



1 + (2fc6 - ko)za 



Inserting these averages in Eq. (|17() and using the explicit 
expression (0} for the correlation function c(fc), we obtain 
for the first cumulant 



2kllk^ 



1 + (2fc(, - kf))zQ 



(ivQ 



2fc2/fco 



1 -I- (2fcb - ko)zo ' 



Here 



(34) 



N 



S{T) = / dk cxp(-rA;)y(z(n))T/(z(n + fc)) 

, (35) 

dk cyij,{-Tk)V{z{n))V{z{n + fc)) 



is the Laplace transform of V {z{n))V {z{n + k)) = s{k) 
and we have used 



dfcc(fc) = — (1-e— ^)"^-^=^. (36) 



Our next task is to calculate s(fc) and S{p). By definition, 
we have 



s(fc) 



j^dl dx dy dz dz' G{x, z; l)5{z - zo)G{z, z'; k)V{z')G{z', y;N -l~k) 



J^dl dx dy dz dz' G{x, z; l)G{z, z'; k)G{z', y; N — I — k) 

I 



(37) 



We adopt the ground-state dominance approximation for This yields for s(fc) 
the overall chain and for the side subchains in Eq. (|37() . 



s(fc)~e'=^« / dzdz'V{z)MzW{z')Mz')Giz,z';k) 
Jo 

(38) 
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and for the Laplace transform 

/>oo 

s{p) ~ / dzdz'v{z)Mz)v{z')M^^')G{^^,^'';p+ kol) 

Jo 

= j dz dz' 5{z- zq)'4>q{z)5{z' - ZQ)ipQ{z')g{z, z' -^p -\- \£q\) 

(39) 



[%l)Q{zQ)Yg{za,zo\p+ |£o|) 

2fc2/fco 



Q{zo,zq;p+ |eo|) 



transform of the Greens function G{z^ z'\ N) with respect 
to the variable N. The function Q{zo, zo;p + \eo\) is cal- 
culated in the appendix. The result is 



G{zo,zq;p) 



sinh (-yBpzo/a) 



(40) 



1 + {2kb - ko)zo' 

where Q{z, z';p) = G{z, z'; N)e~P^ dN is the Laplace Inserting this result in Eqs. H39(l and then (jHU, we obtain 



f3AF, ^ 
N 1 + {2kb - kQ)zo 

X <( /J^A^ Re 



6sinh yerVa^ + fcgzo 



aV6r/a2 + k'^ e\/s^/'^'+'=? - 6/3(^0 + vq) sinh (v/filT^^q^zo) ^ + ~ ''o^^'o 



r 



(41) 



Let us investigate the behavior of this expression near 
the transition point . taking the limit zq — > 0. 
First, we note that in the vicinity of the transition point 
^o' ^ l/^Oi the product kgZQ can be represented as 
kozo = 1 + 5zq. Inserting this expression into Eq. H28|l 
and considering the case zq ^ 1 it is easy to show that 
6 ~ kh- Substituting this representation for fco-^o into 
Eq. |(1T|) . we obtain 



N ~ 



i2kl0^A^ Re 




- 2f3vokb, 

for zo < 1. 
(42) 



D. Determination of the Transition Point 

In the previous sections, we have developed the gen- 
eral theory and obtained expressions for the first cumu- 



lant Ai^i, an approximate expression for the difference 
between the free energies of the original RHP system and 
the reference hompolymer system. Following the idea of 
the reference system approach, we now need to adjust the 
auxiliary parameter wo such that the free energy of the 
reference system approximates that of the original system 
in an optimal way. This implies that at the optimal value 
of uo, the first cumulant Ai^i should be equal to zero: 



(43) 



Looking at the explicit expression for AFi, Eq. 1)41(1 . we 
note that this equation always has one trivial solution 
fcb = or Wo = I'o + Co =0. This corresponds to the 
uninteresting case of a chain which is desorbcd both in 
the original and and in the reference system. The other, 
nontrivial, solution for vq is obtained from the following 
equation: 



J 



/J^A^Re 



6 sinh ^6r/a2 + kl zq 



2kl/ko 1 



aV6r/a2 + k'^ 6^^^/"'+'=' ^« - 6/3(Co + vq) sinh (^6r/a2 + fc^^ zq) ^ + (^^Jj - ^0)^0 F 



pvo = 0. 

(44) 



The transition point for the reference system corresponds nition (|26() of fcg , to 
to the condition zq^o' = ^ (|30|) or, according to the dcfi- 
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Substituting this condition into Eq. yields an equa- 
tion for the transition point ^q*". It will be discussed in 
the next section. 



E. Final Result 

For the sake of convenience, we introduce the dimen- 
sionless variables 
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a 



A = 



/3A 



and 



Our final result for the adsorption transition point 
for heteropolymers (Eqs. (|44|l and H45() ') can then be cast 
into the form 



6 sinh 



(Ver^o) 



$0 -(h =-A'Rc 

/gf e^er £o _ i / j„ ginh ( V6r £q 

(47) 

tr . 

where is the transition point for homopolymers. 

The latter depends strongly on the specific shape of 
the surface potential and in particular on its range zq 

(6/*'^ = 1/620, cf. Eq. However, the shift of the 

transition point for RHPs relative to a homopolymer sys- 
tem, Eq. (|47|) . is much less sensitive to variations of zq. 
In the limit zq <^ 1 the solution takes the elegant form 



^0 -Ch 



(48) 



Figs, n and |21 show adsorption-desorption transition 
curves as a function of the sequence correlation parame- 
ter I/a/F, as calculated from Eq. H47|) for different val- 
ues of Zq. With increasing l/VT, the curves become 
straight lines with the same slope as the asymptotic 
curve, Eq. (|48|l . Only in an interval close to the origin 
does the slope of the curves deviate from the asymptotic 
value. At small zq (zq < 1), the simplified expression H48|l 
provides a good approximation for the full solution (|T7|) . 

The expressions ()47|l and H48|l demonstrate that the 
mean affinity parameter (^0) corresponding to the ad- 
sorption transition is reduced for RHPs, compared to ho- 
mopolymers. In other words, RHPs have a better effec- 
tive affinity to surfaces than homopolymers with the same 
mean affinity. 

Our result also shows that the surface affinity of RHPs 
with exponentially decaying correlations is higher than 
that of RHPs with oscillating correlations. This can be 
seen by rewriting the general expression H48|l for the par- 
ticular cases of purely exponential decay (F e M) 



— tr — 'tr 




and purely oscillating correlations (F G iM.) 



-A^J— ,F iF', F' e 



(49) 



(50) 
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FIG. 1: Shift (?o''-Ch)/A^ of the adsorption-desorption tran- 
sition point for RHPs witli exponentially decaying sequence 
correlations compared to that of homopolymers, vs. correla- 
tion parameter according to Eq. 1471 . The different 
curves show the results for different values of So as indicated. 




l/F 



FIG. 2: 
tions. 



Same as Fig. Q but for oscillating sequence correla- 



III. NUMERICAL APPROACH 

In the previous section, we have developed an analyt- 
ical theory of RHP adsorption based on the reference 
system approach, and subobtained as a main result the 
adsorption-desorption phase diagram, Eqs. l|T7|l and l|I5|l 
and Fig.^ To assess directly the validity of the reference 
system approach, and to test our final result, Eq. H48(l . 
we have also performed numerical calculations based on a 
simple lattice model, following the method of Scheutjens 
and Fleer El. 
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A. Lattice Model 

In the model we use for the numerical calculations, 
polymer chains of N monomers are represented by ran- 
dom walks of length N on the simple cubic lattice. We 
consider polymers which have one end tethered to an im- 
penetrable planar surface. In the z-direction perpendic- 
ular to the plane, the center of mass of a monomer can 
have coordinates z = 0,1, 21, 3/ . . ., where the layer z = 
is directly adjacent to the adsorbing plane. The adsorp- 
tion potential is taken to be short-ranged and sequence- 
dependent and has the following form 

iz,i) _ ( XhU if ^ = ,^ . 

kT " \ if z > ' ' 

i. e., it acts only on the monomers in the layer that is 
adjacent to the adsorbing surface. In Eq. (|51fl Mi denotes 
the type of the i-th monomer in the chain and xm is 
the adsorption interaction parameter for monomers of the 
type M . Note that we consider the simplest case of a 
phantom chain, hence the monomers do not interact with 
each other and the local monomer potential Ui [z, i) does 
not depend on the local polymer concentration. 

One immediately notices one important difference be- 
tween this system and the model considered in the pre- 
vious section: Up to now we have considered free chains, 
whereas we now study tethered chains. However, within 
the ground state dominance approximation used to treat 
chain ends in the last section, the final result for the ad- 
sorption transition will be identical for tethered chains 
and for free chains. Tethering the chains in the numer- 
ical calculations helps to avoid an uncertainty with the 
normalization, because for free chains of finite length one 
would have to introduce a box of finite size. 

We consider two-letter RHP sequences constructed by 
first order Markov chains [Mg = A,B). The underlying 
Markov process is determined by the probabilities to find 
single A and B monomers [J a and fB = ^ — /a^ respec- 
tively), and by the nearest-neighbor transition probabil- 
ities Pi—fj which is the probability that a monomer i is 
followed by a monomer j. It is convenient to introduce 
the correlation parameter 

C = 1 - pA^B - Pb^a (52) 

that characterizes the correlations in the sequence. In the 
case C = one has polymers with uncorrelated monomer 
sequences (BernouUian type), at C > the probability 
to encounter nearest-neighbors of the same type is en- 
hanced, and at C < nearest-neighbor monomers are 
more likely to be of different type (as in alternating 
copolymers). The statistical distribution of the sequences 
is completely determined by the two parameters /a and 
C. 

In the actual calculations, A monomers were al- 
ways considered as adsorbing or sticky {xa is positive), 
whereas monomers of the B type were taken to be ei- 
ther neutral {xb = 0) or repelling from the surface 
{xb = ~Xa)- These two cases will be referred to as 
SN (sticker-neutral) or SR (sticker-repulsive). 



B. Numerical Method 

The statistical weight Gt{z; N) of all conformations of 
tethered chains with one free end in the layer z satisfies 
the following recurrent relation first introduced by Rubin 
|23l | and later used in the more general theory of Scheut- 
jens and Fleer 

Gt{z; iV + 1) = {\Gt{z - 1; iV) + (1 - 2A)Gt(z; N) 
+ \Gt{z + l;N)) 

(53) 



Gt{Q; N+l) = exp(xA/J {(1 - 2\)Gt{Q; N) + AGt(l; N)} , 

(54) 

where A is the probability that a random walk step con- 
nects neighboring layers. On simple cubic lattices, one 
has A = g. Using as the starting point the monomer 
segment distribution G(0; 1) = exp(xMi) and recursively 
applying Eq. H53(l . one can easily calculate Gt{z; N) for 
every chain length N. The statistical weight of all confor- 
mations of tethered chains is then obtained by summing 
over all positions of the free end: 

Nl 

Z = Y,G{z-N) (55) 

The change in the free energy of the tethered chain 
with respect to the free chain in the solution is given by 
/S.F = — log Z. Here the translational entropy of the free 
chain has been disregarded, i. e., the chains are assumed 
to be sufficiently long that it can be neglected. At the 
transition point one has AF = 0, i. e., the energetic ben- 
efit of monomer-surface contacts is equal to the entropic 
penalty connected with the fact that tethering the chain 
at the plane restricts the number of conformations avail- 
able to the chain. 

In order to calculate conformational characteristics 
of the adsorbed chain, one needs to evaluate not only 
Gt{z;n) for arbitrary n — 1,2,...,N, but also a second 
set of functions Gf{z;N — n), which gives the statis- 
tical weight of chain parts between the nth monomer 
and the end monomer N , subject to the constraint that 
the monomer n is fixed in the layer z (whereas the end 
monomer is free). To calculate Gf{z;N — n) one uses 
the same recurrence relations as in Eqs. (|53f) and H54|) . 
but with reverted copolymer sequence and different ini- 
tial conditions G{z; 1) = g -"jv(z)/fcT^ Combining Gt and 
Gf one can then calculate, for example, the average frac- 
tion of A contacts with the surface 

1 ^ 

SA = -^^Gt(0;z)G/(0;A^-i+l)exp(xM.)'5A/..A. (56) 

2 — 1 

Here SMi,A = 1 if Mi = A or otherwise, and the expo- 
nential factor is used to correct for the double contribu- 
tion of the i-th monomer. The normalization constant Af 
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is given by 



JV-l N 



z=0 1=1 



(57) 



The total fractions of A and B contacts can be obtained 
via 



s = -^^AfGt(0;i)G/(0;iV-i + l)e 



Us{0)/kT 



(58) 



In this paper, we present resuhs for chains of length 
N = 500. For every set of model parameters /a and C, 
the adsorption characteristics were calculated as a func- 
tion of the interaction parameter XA for 50 different se- 
quence realizations and then averaged. This corresponds 
to a situation with quenched sequence disorder. In con- 
trast, the method developed by van Lent and Scheutjens 
describes copolymers with annealed disorder. 

A full account of the numerical results will be given 
elsewhere. In the present paper, we are mostly interested 
in the comparison with the analytical theory, and in par- 
ticular, in the two questions that are discussed in the next 
section. 



IV. COMPARISON OF THE TWO 
APPROACHES 

A. Is the reference system approach vahd? 

Introducing the reference homopolymer system in Sec- 
tion IIIBI we have assumed that its conformational and 
thermodynamic properties are reasonably close to those 
of the original RHP system. Later, the optimal value of 
variational parameter vq has been chosen according to 
the requirement that it minimizes the free energy differ- 
ence between the original and the reference systems ()43f) . 
Therefore, the free energies of the two systems are auto- 
matically close. However, one cannot be sure that this 
choice of va guarantees good correspondence between the 
conformational characteristics. This shall be tested first. 

One of the most important and indicative conforma- 
tional property of the adsorbed chain is undoubtedly the 
total fraction s of adsorbed monomer units. It is plot- 
ted in Fig. 121 for RHPs with different composition and 
structure, and for homopolymers with the same free en- 
ergy. One can see that near the adsorption transition 
{AF = 0), the value of s for RHPs is close to that for 
homopolymers. The "worst" RHPs are those which con- 
tain a majority of surface-repelling monomers. In the 
strong adsorption regime, however, the values of s differ 
considerably. 

This observation can be explained by the fact that in 
the vicinity of the transition point, chains have relatively 
few surface contacts and long loops and tails. Such con- 
formations can be realized equally easily by both ho- 
mopolymers and RHPs, in spite of the "gaps" of neu- 
tral or repelling monomers in the RHPs. Since only few 
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FIG. 3: Fraction s of adsorbed segments for homopolymer 
and RHPs of different type (SN or SR) and composition (/a) 
as a function of the free energy. For each of RHP type and 
composition curves for C = 0, 0.25, 0.5 and 0.75 are shown. 



units adsorb to the surface, the system can avoid en- 
tropically and energetically unfavorable contacts of neu- 
tral/repelling units with the surface. In the strong ad- 
sorption regime, it is thcrmodynamically more favorable 
for homopolymers to have many contacts with the sur- 
face. This is not necessarily the case for RHPs, because 
non-adsorbing monomer units still avoid to get involved 
in the interactions with the surface. 

In the strong adsorption regime, the reference system 
approach of section lTlBl thus becomes questionable. How- 
ever, it seems to be applicable in the vicinity of the tran- 
sition, and should therefore be suitable to study the tran- 
sition point. 



B. Is the analytically predicted phase diagram 
correct? 



Next we check the scaling relation, Eq. H48|) . describ- 
ing the shift of the adsorption transition point of RHPs 
relative to homopolymers. 

To compare Eq. H48|) with the numerical results, we 
need to do two things: (1) find the transition points for 
lattice RHPs with different composition and "structure" 
and (2) "translate" these results into of parameters of the 
continuum model parameter {^q, A, and F). 

The transition point Xa lattice model was deter- 

mined numerically from the condition AF = for every 
single realization of a RHP. Then the average over the 
values x^^ was performed. 

The next step is to find the relation between the sta- 
tistical parameters of the continuous model and the dis- 
crete lattice model. For homopolymer systems this prob- 
lem has been solved by Gorbunov et al. According 
to this work, the coefficient connecting the continuous 
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and the discrete polymerization degree is unity, and the 
length a in the continuum model corresponds to the lat- 
tice spacing I in the case of a simple cubic lattice. To 
find the mapping for ^o, and F, we calculate the mean 
monomer-surface interaction parameter 



Xo = IaXa + (1 - fA)XB 



(59) 



and the monomer-monomer correlation function in terms 
of the monomer interaction energy. The calculation is 
straightforward and yields 



Cx W ={ix{x) ~ Xo) {x{x + n)- xq)) 



(60) 



-0.4 




If we represent the correlation function in an exponential 
form like Q 



A 7r 

7. 



Cx(ri) = A^exp(-r^n) , 
we obtain the equivalent parameters A and T 

r^ = -iog(c) 



A, 



XB 



)VfA{l-fA). 



(61) 



(62) 



Since the units for the polymerization degree are the same 
in both models, this yields F = F-^. 

The two remaining parameters and A are related to 
the monomer-surface affinity. The mapping of this "en- 
ergy scale" is much less straightforward. We shall adopt 
a procedure suggested in Ref. [3 and adjust the energy 
scale such that the slope of the monomer density pro- 
files of adsorbed homopolymers at the surface is the same 
in the vicinity of the adsorption transition. Comparing 
Eq. (15) of Ref. 13 with the zq ^ hmit of Eq. we 
obtain 



XH - Xh 6 



(63) 



for homopolymers on simple cubic lattices. We assume 
that the factor (5/6) thus relates the energy scales in our 
two models, i. e., the continuum model and the lattice 
model. Thus we conclude 



(rtr ttr 



Xo 



A 

a: 



(64) 



The approach has the drawback of being rather indirect. 
Moreover, the energy mapping applies, if at all, only in 
the close vicinity of the transition point. 

Inserting these relations for T- and A into Eq. H48|l 
we obtain the following relation between equivalent lat- 
tice parameters 



xr- 



tr 

Xh 



(65) 



FIG. 4: Shift Xo" — Xh oi the adsorption transition in RHP 
systems relative to homopolymer systems for the lattice model 
vs. A^/r^''^. The solid line shows the best linear fit. 



Thus Xo" ~ X*^ is expected to depend linearly on A^ /V^, 
with a slope -5/^6 w -2.04. 

The results for X([ ~ Xh are shown in Fig. 0] One 
can see that the numerical calculations reproduce well 
the expected linear dependence between X([ ~ X% and 
A^/Vr. However, the slope (w —0.754) is different from 
that predicted by Eq. 



V. CONCLUSIONS 

To summarize, we have studied the adsorption of sin- 
gle ideal heteropolymer chains onto homogeneous planar 
surfaces. After applying the replica trick we have intro- 
duced a reference homopolymer system chosen such that 
the thermodynamic properties were as close as possible 
to those of the original RHP system. This approach al- 
lowed us to treat the problem analytically and to obtain 
the adsorption-desorption phase diagram of the RHP. In 
particular, we have considered the case of RHPs with 
primary sequence correlations described by a de-Gennes 
correlator. We found that the adsorption is enhanced 
with increasing strength of correlation (A) and correla- 
tion length (1/F), both for the cases of decaying and os- 
cillating correlations. Our result could be summarized 
in a very simple equation for the shift of the adsorption 
transition in correlated RHPs relative to homopolymers, 
Eq. gHIl- 

To test the reference system approach, numerical lat- 
tice calculations were performed. These results show 
that in the vicinity of the transition point, our varia- 
tional is well justified. Furthermore, the calculation of 
the adsorption-desorption phase diagram confirmed the 
analytical scaling prediction. However, we also found 
that far from the transition point, one cannot "construct" 
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a reference homopolymer system that is good in both 
a conformational and a thermodynamical sense. A ref- 
erence system perhaps more adequate to describe this 
case could be a multiblock-copolymeT system with block 
lengths equal to the average length of A and B blocks in 
the random copolymers. In the lattice model, these are 
equal to sa = fA)(l - c)] and sb = -c)], 

respectively [23 ). This could be interesting for future 
work. 

We found that our main result, Eq. (|48|l . describes 
the numerical result very satisfactorily from a qualitative 
point of view. Unfortunately, true quantitative agree- 
ment could not yet be established. The reason lies partly 
in the fact that the identification of energy scales in the 
lattice model and in the continuum model is not self ev- 
ident. We have tested a mapping procedure which ad- 
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justs the monomer profiles for adsorbed homopolymcrs. 
Within this approach, the theory seems to overestimate 
the shift of the adsorption transition in RHP systems. 
However, the mapping can be questioned, as we have dis- 
cussed in the last section. Furthermore, we are compar- 
ing a theory for RHPs with continuous (Gaussian dis- 
tributed) monomers with numerical calculations for two- 
letter RHPs, which might also lead to problems. 



Acknowlegdement 

The financial support of the Deutsche Forschungsge- 
meinschaft (SFB 613) is gratefully acknowledged. 



Appendix: Calculation of G(zo, zo;p) 



Our starting point is Eq. H18|l with the initial condition 119|l and the boundary conditions H2U|) . Laplace transforming 
with respect to N gives 



-pg{z,z';p) + G{z,z';0) 



6 dz^ 



Pwo5{z - zo)g{z, z';p), 



(66) 



where g{z, z';p) is the Laplace transform of G{z, z'; N) with respect to N, g{z, z';p) = G{z, z'; N)e P^dN. Taking 
into account the initial condition (|19|1 . we obtain 



-pg{z, z';p) + 5{z - z ) 



a? d^g 

6 dz"^ 



■ Pwo5{z - zo)g{z, z';p) 



with the boundary conditions 

g{z,z';p)\^^oc =0, g{z,z';p)\^^o = 
Laplace transforming H67|l with respect to z with H2l)|) gives 



(67) 
(68) 



-pgis, z'-p) + e-^'« = -- [s^g{s, z'-p) - s5(0, z';p) - ^'(0, z'-p)] - ^w^gizo, z';p)e-'°' 

where g{s,z';p) is the Laplace transform of g{z,z';p) with respect to z, g{s,z';p) ~ g{z, z';p)e~''^dz, and g' 
stands for dg/dz. With the boundary condition at z = H68|) we obtain 



-pg{s,z';p) + e -- 
Solving this equation for g(s, z';p) yields 



[s^g(s, z';p) - g'{0, z';p)] ~ f3wng{zo, z';p)e~ 



g{s,z';p) = 



g'iO,z';p) - (6/3u;o/a')e(zo,^';p)e-"°^ - i6/a^)e- 



— {(!>p/o?) 

Taking the inverse Laplace transform of (|70|l with respect to s gives 

6/3iyo . , 
— - — smh 



e(z,z';p) - -^{sinh(v^^^) g'{Q,z'-p) 
6 



■ sinh 



Qp - 



where 9{z) is the Heaviside function 



1, if z > 
0, if z < 



(69) 



(70) 



6p - 



g{zo,z';p)e{z - zo) 



(71) 
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Setting z ~ zq in Eq. (|71|l . we obtain 



g{zo,z';p)^ jsinh (^^6^^) ^^'(0, z'; p) - ^ sinh ^^Gp "^^-^^ d{zo-z') 



(72) 



and then insert G^zq, z']p) back into Eq. (|71|) (second term): 
c;(z,z';p) = -^|sinh(v/6^^) C/'(0,z';p) 



sinh v6p 



sinh I 



(v^f) t;'(0,z';p) 



^ sinh ( -^Z 6p 



a 

z — z 
a 



■ sinh -v/ 6p 



0{z-zo) 

Zq - z') 



(73) 



To find the unknown C7'(0, 2';p), we require Q{z — > oo,z';p) to vanish and z' to be finite. This imphcs that the 
coefficient of the growing term term e'^'-l'^ must be zero. Using (|73|l . this requirement leads to the following result 
for e'(0,z';p): 



g e-^^^'/o _ Mho sinh |^^£o^^ e-^^o/"6'(zo - z') 
^ (0' ^ ;P) = ^ 1 _ Mho gjjjl^ ( V6p f ) e- V6F-o/'^ 

To find t/(zo, zo;p), we substitute H74|l into H73|) and set z = zq and z' = zq. The result is 



(74) 



^(zo,zo;p) = -|= 



sinh (^/6pzo/a) 



a2 ev^-o/a _ M^_a^ sinh (V6^zo/a) 
sinh (y^zo/a) 



(75) 



ai/p/eev^^o/a _ sinh (^^zo/a) 
I 
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